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ABSTRACT 

Stellar population studies show that low mass galaxies in all environments exhibit 
stellar halos that are older and more spherically distributed than the main body of the 
galaxy. In some cases, there is a significant intermediate age component that extends 
beyond the young disk. We examine a suite of Smoothed Particle Hydrodynamic 
(SPH) simulations and find that elevated early star formation activity combined with 
supernova feedback can produce an extended stellar distribution that resembles these 
halos for model galaxies ranging from W200 — 15 km s~^ to 35 km s~^, without the 
need for accretion of subhalos. 
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1 INTRODUCTION 

The metal-poor stella r halo is th e oldest extended struc- 
ture in the Galaxy (|Oortl Il924l ). Its internal structure, 
kinematics, and metallicity encode information about the 
early assembly of the Galaxy. Inspired by the constraints 
available in the Milky Way, several groups have be- 
gun intensive observational studies of comparable struc- 
tures in other large spiral galaxies, most notably in 
M31 and M33 (iFerguson et al.l I2OO2I: iGuhathakurta et al.l 



the accret ed stars (cf., Bekki et al.| 20031: Brook et al.l 
2004 Bu Uock fc JohnstonI l2005l: Renda et all l2005l: 
iBckki fc Chibal l2005l: jGauthier erahl l2006l : lAbadi et al.] 



I2OO5I : [Gilbert et al.ll2006l : llbata et al.ll2007l ). in several ded 
icated HST programs (e.g., see p reliminary results in 
ISeth et al.l l 2007l : lde Jong et all l2007l), and in isolated fields 
in nearby galaxies llHarris fc Harris! I2OOOI: iMouhcine et al.l 



l2005l : lTikhonovl2006al : ISeth et al.H2007l : lde Jong et al.ll2007l ). 

These studies use resolved stellar populations to trace the 
extended structure of galaxies down to surface brightnesses 
that are an order of magnitude fainter than what is possible 
with ground-based surface photometry. 

Both observations and simulations suggest that the 
stellar halos of massi ve galaxies are a rem nant of early 
hierarchical merging (|Searle fc ZinnI 1 19781 ). As galaxies 
grow, they repeatedly accrete smaller galaxies. Dynamical 
friction pulls the accreted galaxies inwards, while tidal 
disruption strips the galaxies' stars. These stars are left 
on orbits with large apocenters, creating an extended 
stellar halo. Increasingly mature numerical simulations 
now make specific predictions for the profile of stellar 
halos, their degree of substructure, and the metallicity of 
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. -_^_Jg , 

I2OO6; Font et aLlbOOQ : ISckki 2008^ 

Less theoretic al study ( cf., iMashchenko et ahl I2OO8I : 

lBekki,.2008 : Valcke et al.ll2008l ) has been devoted to the ex- 
tended stellar structure of dwarf galaxies, in spite of the fact 
that old stellar populations have been observed surrounding 
every low mass galaxy that has been studied to sufficient 
depth. We present a review of these observations in Sj^j At 
first glance, the observations are consistent with dwarf halos 
being analogous to the halos seen in the MW and other mas- 
sive spirals, suggesting that stellar halos may be a common 
by-product of galaxy assembly at all mass scales. 

In spite of their similarities, significant differences be- 
tween dwarf stellar halos and massive galaxy halos are be- 
coming apparent, suggesting that they may not be true 
analogs of one anotherQ First, deep color-magnitude di- 
agrams (CMDs) constructed just outside t he main opti- 
cal bo dy of the Phoenix "tra nsitional" dSph dHidalgo et abl 
l2003j ) and the Leo A dirr (|Vansevicius et all |2004| ) indi- 
cate a non-negligible population of 2-8 Gyr old stars at these 



^ In spite of these differences, we will continue to refer to the 
low surface brightness, extended distribution of older stars sur- 
rounding dwarf galaxies as a "halo" since it best describes the 
morphology of the extended stellar population. However, the kine- 
matics, metallicities, age distributions, and origins of the stars in 
the dwarf "halo" are not necessarily exact analogs to the halo of 
the Milky Way. 
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larger radiiEl This "intermediate age" population has no 
counterpart in the ancient halo of the Milky Way or the 
outer halo of M31 (jCuhathakurta et all [20051 : iKalirai et all 
[2006 ). although such a population does exist in the inner 
halo of M31, which is l i kely c ontaminated by thick disk 
stars ([Brown et alj [200^ . l2008h . Second, the surface den- 
sity of dwarf galaxies' stellar halos seem to fall off expo- 
nentially, rather t han as a power-law (see radial profiles in 
Vansevicius et a l] [2004 : Hidalgo et al.|[2003d : lAparicio et alJ 



19971 . for example) . The somewhat brighter low mass galaxy 
NGC 4244 (Ms = -17.0) also appears to host a radially ex- 
ponential distr ibution of halo s tars, although a power law is 
not ruled out (|Seth et al.ll2007l '). 

If we apply standard hierarchical structure formation 
arguments to form dwarf stellar halos, it is unclear how 
these features would be produced. First, it is unlikely that 
any of the low mass sub halos that dwarf galaxies accreted 
contained sufficient stellar material to populate the halos 
jPurcell et al.ll2007l ). Second, since most merging happens 
early in the history of the Universe, it is unlikely that merg- 
ing can routinely produce significant intermediate age stellar 
populations. 

To understand how the mass-dependent differences in 
the properties of stellar halos might arise, one should con- 
sider alternate possible origins for the dwarf galaxy halos. 
In this paper, we explore possible mechanisms for halo for- 
mation focusing on those that are intrinsic to the galaxies 
themselves rather than dependent on mergers or interac- 
tions. 

For isolated galaxies, stars located in an older stellar 
halo could have formed m situ in a contracting envelope of 
star formation or could have formed in a central region and 
then been ejected. Using our N-body simulations, we demon- 
strate that these mechanisms can indeed produce extended 
stellar halos in isolated dwarfs, without the need to accrete 
stars from merging subhalos. We show that the resulting ha- 
los have exponential profiles and age gradients whose char- 
acteristic scale lengths agree with observations. 

This paper is structured as follows. In fj^l we present 
observational evidence for dwarf halos. In fjS] we describe 
the model we use for dwarf galaxies that form stellar halos. 
In 33] we examine the age gradients, stellar density profiles, 
star formation histories, and halo shapes that are produced 
in the models. 



2 EVIDENCE FOR DWARF GALAXY HALOS 

Over the past decade, access to wide-field optical and near- 
infrared imaging has increased dramatically. These techno- 
logical advances have expanded the opportunities to obtain 
wide-field spatially-resolved CMDs of the stars within dwarf 
galaxies. Access to similar imaging capabilities in space, 
through the Hubble Space Telescope, has allowed compa- 
rably high-quality CMDs to be derived well beyond the con- 
fines of the Local Group, out to ~5 Mpc. 

This growing body of observations has led to wide- 
spread confirmation of what had already been suspected 

^ Such an intermediate age population could not yet have been 
detected in other dwarf galaxy CMDs, which are typically much 
shallower than the observations of Phoenix and Leo A. 



from ground-based surface photometry, namely, that star- 
forming dwarf galaxies show population gradients in their 
outer regions. The CMDs of dwarf irregular and "transi- 
tion" dwarf galaxies universally show extended envelopes 
dominated by older red giant branch (RGB) stars. These 
envelopes can extend many times further than the young 
bright blue stars, and often have different axial ratios and 
morphologies than the embedded young component or the 
extended distribution of HI (e.g., NGC 6822, IC 1613; 
iBattinelh et al.ll2006l : [Bernard et al.ll2007l '). 

A thorough review of the literature yields more than 
50 papers on ~ 40 star forming dwarf gala xies hosting an 
extend ed halo of red gian t bran ch stars ([van Dvk et al.l 



Alonso-Garcfa et aL , 

200(i 



20061 : [Aparicio et 



Aparicio et al.l 



In" 



1997 



■ .200C 
Battinelli fc De mersI l2004allb 



Demers et al.l 



2004 



119981 
|2006| 



Drozdovskv fc Tikhonov[ [2000l: [Callart et al 



Aparicio fc Tikhonov 
Battaglia et al.' '2006": 

2005, 2006; Battinglli ct al. 2006, 200 7.: [Bernard et al.[ 
20071: ICaldwell et all 



[de Jong et al 



iDoublier et al.l 



2008 



200C 



2004 



Han et all [ 19971: [H eld et al' '200l'. '1999*: 'Hidalgo et al 



Lee & Bvun|[l999l: Lee fc Hwang 2005: 


Letarte et al.l [20021: 


Martfnez-Delgado et al.l 1999i. ,999ai: 


McConnachie et al. 



2009; iMiller et aL. 
200i; [Minniti fc Ziilst: 
Noel fc Gallart .2007. : .Piersimoni e t al 
2000 !; 



20011: iMinnitil I2OOII: iMimiiti et~aL 
ilstra' '1997'; 'Minniti e t al.[ [l999l: 



1999 : Reikuba et alj 



Savia ne et al. 2000l : Schu lte-Ladb eck et al 



199E 



20021: iTikhonov fc Galazutdinova' '200 5: ITikhonov[ [2006a[ lb 



Vansevicius et a l. 2004; Ziilstra fc Minniti 19991 ). Earlier 



re views can be fou n d in iSaviane et al.l (|200ll ) and in part 



of [Kunth fc 6stlin[ ([2OOOD . There were no cases reported 
in the literature where this effect was looked for, but 
not found[3 In Figure [T] we show the properties of the 
galaxies with extended RGB distributions as a function 
of a bsolute B magnitud e and tidal index, both taken 
from iKarachentsev et al.l (|2004l )'s Catalog of Neighboring 
Galaxies. The tidal index indicates the strength of ex- 
pected tidal forces on the galaxy, based on the nearest 
massive neighbors; negative values indicate galaxies that 
are essentially isolated. 

Figure [1] shows that the presence of an extended halo of 
older stars is seen in all environments over a range of stellar 
masses. The literature thus suggests that the formation pro- 
cess must be intrinsic to the galaxies themselves rather than 
due exclusively to tidal effects from larger galaxies. Unfor- 
tunately, no data yet exist on the fraction of stellar mass 
in the extended halo, and thus we can say only that such 
envelopes exist, and are ubiquitous in low mass dwarfs. 

Stellar population gradients also are found in the 
majority of dwarf spheroidals with sufficiently ex- 
tended observations (Sculptor, LeoII, Audi, Andll, 
NGC 185, NGC 147 ESQ 4 10-G005; [Tolstoy et al] 
'2004 1: iKomivama et al] [20071: 
McConnachie et al. 2003; 



Da Costa et al.[ 

Martmez-Delgado et al.l 



Mcuonnacme et ai.l l^^UUYl: iMartmez-ueigado et ai.l iitJija 
Hanet al.l [ 19971 : [ Karachentsev et al.[ [20001 ). Theseobserva- 



^ The one reported non- detection in NGC 31 09 bv lDemers et all 
1I2OO3I ) is contradicted bv lMinniti et all l ll999l) . who found an ex- 
tended RGB envelope using data that reaches a full magnitude 
deeper. 
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Figure 1. Absolute _B-band magnitude and tidal index for dwarf 
galaxies with repor t ed RG B halos. Quantities are taken from 
[Karachc ntsev et al. Catalog of Neighboring Galaxies. 

The tidal index indicates the degree of expected tidal influence 
from nearby massive galaxies, with negative values indicating 
truly isolated systems. Open circles are Local Group galaxies with 
the Milky Way or M31 as the primary perturber. Solid squares 
are galaxies in the M81 group. Solid circles are all other galaxies. 
For reference, the asterisk in the upper right corner indicates the 
properties of the Small Magellanic Cloud. 



tions are more subtle to interpret, due to the lack of current 
star formation. Instead, the population gradients typically 
manifest as changes in horizontal branch morphology with 
radius, such that bluer, more extended horizontal branches 
(favoring older ages and/or lower metallicities) are seen 
towards the outskirts. The gradients can also be traced 
using intermediate age asymptotic giant branch (AGB) 
stars, which are found to be more centrally concentrated. 
Unlike the dwarf irregulars and transition galaxies, however, 
there are 4 well-studied dwarf spheroidals (Leol, Carina, 
Tucana, and possi bly Draco) in which gradients have 
not been detected llHeld et all | _2000!; 'S aviane et all 1 19961 : 
iDa Costa et al.l I2OO0I : ISegaU et a l. 2007). However, given 
that the Ukely formation scenarios for dwarf spheroidals 
involve extensive dynamical threshing o f dwarf irregulars a s 
they are accreted by larger halos (e.g.. iMaver et al]|2007l ). 
it is more surprising that any dwarf spheroidals could have 
retained population gradients imprinted before they fell in. 

While the evidence from color-magnitude diagrams 
confirms the presence of older stars in the envelopes of 
dwarf galaxies, few observations yet constrain exactly how 
old those stars are. A simple stellar population forms a 
prominent red giant branch after only a Gigayear or so, 
and thus some of the stars in the RGB envelopes might 
be of intermediate age, provided that the main sequence 
turnoff is too faint to be detected in the existing obser- 
vations. Observations of galaxies beyond the Local Group 



rarely reach the horizontal branch, and therefore are in- 
capable of detecting main sequence turnoffs of younger 
than ~ IGyr. The RGB envelopes may therefore host 
an intermediate age stellar population, and may not nec- 
essarily be strict analogs of the ancient stellar halo of 
the Milky Way. Indeed, extensive near-infrared observa- 
tions of AGB stars in the halos of dwarf galaxies demon- 
strates the presence of intermed iate age (1-lOG yr) popu- 
lations f e.g.. iDemers et 'al.,.2004: .Battinelli fc Demcrs 2006I. 
l2004al lbl: lOemers et al.1 12OO6I : iLetarte et al.ll2002l ). Such in- 
termediate age populations have also been verified by direct 
detection o f main sequence turno ffs in the transition dwarf 
Phoenix bv lHidalgo et all l|2003bl ). although the fraction of 
intermediate age stars is small. Given the limited evidence 
to date, one should not assume that the "halos" of red giant 
branch stars around dwarfs are purely ancient. 

Taken together, the extensive body of literature above 
indicates that extended envelopes of red giant branch stars 
and strong radial population gradients are common in 
dwarf galaxies, with the exception of some subset of dwarf 
spheroidals in the Local Group. The gradients indicate that 
the outer envelopes are older than the central galaxy, but 
that they may contain some fraction of intermediate age 
stars. The ubiquity of the phenomena in galaxies of all envi- 
ronments indicates that halos can be produced by phenom- 
ena internal to the galaxy, rather than being induced by 
tidal interactions with massive companions. The presence of 
intermediate age stars in some halos further indicates that 
dwarf halos are not necessarily a by-product of processes 
that operate only during the initial formation of the galaxy. 



3 SIMULATIONS 

To examine dwarf halo formatio n, we begin using the same 
suite of simulations described in IStinson et al.l l|2007l ). The 
model galaxies vary in total mass from lO'^M© to 1O^°M0 
and have velocities at the virial radius (i'20o) of 15, 20, 
25, and 30 km s~^, resulting in stellar velocity dispersions 
comparable to Local Group dwarfs. Each model galaxy ini- 
tially consists of a sphere of 10^ gas particles in hydrostatic 
equilibrium with a sph ere of 10^ dark matter particles, ini- 
tially distributed with a iNavarro et al.l (|l997l ) density profile 
{c — 8). The baryons account for 10% of the total galaxy 
mass. Dark matter particles in ev ery model are left with 
their random equilibrium velocities (|Kazantzidis et al ]|2004l ) 

while the gas is spun with a uniform circular velocity so that 
1 

— 0.039, where jgas is the average specific an- 



X = 



-3- 



gular momentum of the gas, and E and M are the total en- 
ergy and mass of the halo. Although these initial conditions 
are not ideal analogs of the likely mechanisms for gas accre- 
tion, the mechanisms we explore in this paper only manifest 
after the formation of a settled star forming system. Thus, 
the properties of the resulting halos are unlikely to depend 
sensitively on transient efi'ects linked to the specific choice 
of initial conditions. 

We supplement these models with a higher resolution 
simulation of the 15 km s~^ galaxy to maximize the number 
of particles in the diffuse stellar halo. The high resolution 
model is the same in every respect to those described above 
except that it contains ten times more gas and dark matter 



4 Stinson et al. 



particles. Each gas particle in the high resolution simulation 
initially contains 140 M© of material, and each star particle 
represents an initial mass in stars of 40 M©. 

To these, we add simulations with varying baryon frac- 
tion. In a cosmological context, reionization has a strong 
impact on small galaxies and puts an effective mass limit 
on galaxies that can retain enough gas to form stars 



(Quinn ct al. 1996: Ha ardt fc Madaulll996l : lBarkana fc Loebl 
11999 ; Bullock et al. 200oh . We have not included reioniza- 
tion explicitly, and instead capture the relevant physics 
by simply excluding increasing fractions of gas from the 
galaxies as would be expected if a UV background pre- 
vented some fraction of the gas from cooling into a disk. 
A more realistic simulation that included the full cosmo- 
logical context would be required to verify this scenario, 
but resolving s tructure at this sca l e wou ld be prohibitively 
expensive (cf. iM ashchenko et al] (|200 8l) require 10® pro- 
cessor hours to evolve such a dwarf to z — 5). Empir- 
ically, low mass galaxies have high mass to light ratios, 

and at least an order of magnitude scatter in their inferred 

' II 1 I II 1 

baryon fr action (van den Bosch ot al. 2003; Lin ct al. 2004!; 

Yang et ^ l. 2005: Hoekst r a et al.. .2005. : .Mandelbaum et al.i 
20061 : IZaritskv et aDl2007l l. To mimic this, we simulate three 
35 km s~^ galaxies with baryon fractions (/[,) of 10 %, 
1%, and 0.5%. The varying galaxies contain 10^ gas 
and 10^ dark matter particles, but all have a total mass o f 
1.365 X 1O^°M0, the mass above which iQuinn et all (119961 ) 
find a galaxy can survive reionization («200=35 km s~^). 

The si mulations were run using the parallel SPH code 
GASOLINE (iWadslev et ahlliooih . GASOLINE solves the equa- 
tions of hydrodynamics, and includes radiative cooling. The 
cooling assumes ionization equilibrium, an ideal gas with 
primordial composition (metal line and H2 cooling are not 
included), and solves for the equilibrium abundances of 
each ion species. The scheme us e s the coUisional ioniza- 
tion rates reported in lAbel et all lll99'i1) . the radiative re- 
combi nation rates from lBlackl 1 198lh and lVerner fc Ferlaridl 
(|l996t ). bremsstrahlung, and line cooling from ICenI l| 19921 ). 
The energy integration uses a semi-implicit stiff integrator 
independently for each particle with the compressive heat- 
ing and density (i.e. terms dependent on other particles) 
assumed to be constant over the timestep. Gravity is cal- 
culated for each particle using tree elements that span at 
most 6 — 0.7 of the size of the tree element's distance from 
the particle. GASOLINE is multistepping so that every par- 
ticle's time-step Atg^av = ^\f^^ where 77 = 0.175, is 
the gravitational softening length, and is the accelera- 
tion. For gas particles, the time-step must also be less than 

Atgas = J^Courant JYqr^J^-^f^^T^^J^, wherC 77Courant = 0.4, hi is 

the gas smoothing length, a is the sound speed, a = 1 is 
the shear coefficient, /3 = 2 is the viscosity coefficient, and 
fJ-MAx is described in full detail in'Wadsl ev et al] (|2004l ) . 

The star formation and feedbac k recipes are the " blast- 
wave model" described in detail in lStinson et al.l l|2006l ). and 
they are summarized as follows. Gas particles must be dense 
(wniin = O.lcm""^) and cool (Tmax = 15,000 K) to form stars. 
A subset of the particles that pass these criteria are ran- 
domly selected to form stars based on the commonly used 
star formation equation. 



dt 



tdyn 



(1) 



where M* is mass of stars created, c* is a constant star for- 
mation efficiency factor, Mgaa is the mass of gas creating the 
star, dt is how often star formation is calculated (1 Myr in 
all of the simulations described in this paper) and tdyn is the 
gas dynamical time. The constant parameter, c* is tuned to 
0.05 so that the simu lated Isolate d Model Mi l ky Wa y used in 
IStinson et al. I (j2006) matches the lKennicuttI l| 19981 ) Schmidt 
Law, and then c* is left fixed for all subsequent applications 
of the code. This star formation and feedback treatment w as 
one of the keys to the success of iGovernato et al] l|2007h in 
producing realistic spi ral galaxies in a cosm ological simula- 
tion and the success of lBrooks et al.l (|2007h in matching the 
observed mass-metallicity relationship. 

At the resolution of these simulations, each star particle 
represents a large group of stars (40-3500 M0). Thus, each 
particle has its stars partitioned int o mass bins based o n 
the initial mass function presented in iKroupa et all l|l993l ). 
T hese masses ar e correlated to stellar lifetimes as described 
m iRaiteri et all (|l996h . Stars larger than 8 Mq explode 
as supernovae during the timestep that overlaps their stel- 
lar lifetime after their birth time. The explosion of these 
stars is t reated using the analytic m odel for blastwaves pre- 
sented in McKee fc Ostrikerl (Il977l ) as described in detail in 
IStinson et all l|2006l ). While the blast radius is calculated 
using the full energy output of the supernova, less than 
half of that energy is transferred to the surrounding ISM, 
EsN = 4 X 10^° ergs. The rest of the supernova energy is 
radiated away. 



4 RESULTS 

In this section, we examine the spinning gas halos that cool 
radiatively into centrifugally and pressure supported star 
forming disks. First, we provide evidence for their extended 
exponential stellar density profile and the presence of an age 
gradient in the models. We then examine the mechanisms 
responsible for these features. 



4.1 Density Profiles 

In Figure [2] we show the radial surface density profile of the 
stars formed in the simulations for several difi'erent mass 
simulations. The halo surface density profiles generally fol- 
low an exponential, though there are breaks where the scale 
lengths change. We fit these profiles with a double expo- 
nential. The inner and outer exponential scale lengths are 
reported in Table [T] along with the break radius that sep- 
arates each exponential segment in Figure [2] Similar bro- 
ken exponential profile s are evident i n observ ations; for ex- 
a mple, see Figure 1 in Hidalgo et al.l l|2003bl ) and Figure 3 
in IVansevicius et al.l l|2004h . Alternatively, the surface den- 
sity profiles can be fit with a series of power laws, with a 
flat inner core, a middle region that declines as r~^, and 
an outer region that declines steeply as or steeper. Al- 
though the two component exponentials are reasonable fits 
to the overall profile, there are significant deviations that 
are more pronounced in the lower mass galaxies. If mecha- 
nisms involving supernova feedback create these outer halos, 
it makes sense that deviations would be more pronounced in 
shallower potential wells where the supernovae have a larger 
relative effect on the gas dynamics. Figure [2] also shows 
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Figure 2. Projected edge-on radial stellar density profiles of the 
simulated galaxies. Stars are averaged over radial bins that each 
contain equal numbers of star particles. The lines correspond to 
the best two piece exponential fits to each radial set of points. 
The more massive halos have higher central surface densities and 
longer exponential scale lengths. Note that these plots do not 
show the number of stars brighter than some magnitude limit (i.e. 
Figure |9]l , and thus they are not exact analogs of observational 
studies. 




0.0001 



2 4 
R[kpcl 



Figure 3. Projected face-on and edge-on radial stellar density 
profiles of the high resolution 15 km s~^model. Stars are averaged 
over radial bins that each contain equal numbers of star particles. 
The lines correspond to the best quadruple piece exponential fit 
to each radial set of points. In the face on profile, the segments 
ending at 0.49, 0.8, 1.22, 2.22, and 4.5 kpc have scale lengths of 
0.35, 0.52, 0.30, 0.49, and 0.81 kpc, respectively. The edge-on view 
segments ending at 0.49, 0.8, 1.22, 2.08, and 4.5 kpc have scale 
lengths of 0.29, 0.41, 0.27, 0.49, and 0.82 kpc. 



that the stellar scale l enRth increases with mass, as expected 
from analytic models l|Dalcanton et al.|[l997l : lMo et al.ll 19981 : 
Ivan den Boschll200ll ;). 

Figure [3] shows face-on and edge-on projections of the 
surface density profile of the high resolution 15 km s ~^ 
model. The high resolution profile can be modeled with a 
series of exponential segments with scale lengths ranging 
from 0.3 to 0.8 kpc for both projections. The outer scale 
lengths for the two orientations become more similar at large 
radii suggesting an increasingly isotropic, spherical shape 
with radius. We return to this point in i|4.4l 



4.2 Age Gradients 

CMDs of dwarf galaxy halos show well populated red giant 
branches with little evidence for bright main sequence stars 
(32]). These data suggest that the halos are systematically 
older than the central star forming regions. We identify a 
similar age gradient by examining the distribution of stars 
in our models three different ways. We first show the mean 
age gradient (Figure |4]). We then look in detail at the full 
distribution of stellar ages as a function of position (Figures 
[Slini[71 and [5}. Finally, we generate model CMDs from our 
simulations within radial bins (Figure |9}. 

Figure |4] shows the mean age profile of stars as a func- 
tion of radius. In all the halos, stars at the smallest radii have 
the youngest mean ages. However, the mean age in the center 
of the galaxies is never less than 6 Gyr, the mean age ex- 
pected if stars have been forming continuously for 12 Gyr. In 
the lowest stellar mass halos (15 km s"^ and /6=0.005), the 
stars at large radii have the oldest mean ages. In high stel- 
lar mass halos, the age profile turns over to slightly younger 
ages at larger radii, as might be expected for classic inside- 
out disk growth. Lower mass galaxies do not show this trend 
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Figure 4. The mean stellar age versus radius. The mean is taken 
of stars grouped together into 40 radial bins with equal numbers 
of particles. The thin lines represent the models run with a baryon 
fraction of 10% and variable mass. The thick lines represent the 
models that have a constant mass corresponding to v^i,- = 35 km 
s~^ and variable baryon fraction. 



because they lack a stable star forming disk. In all cases, 
however, the outer regions are quite old (> 8 Gyr). 

Figures [5] and [6] provide a more detailed examination 
of the behavior seen in Figure [4] The plots show the final 
location of stars versus their formation time. They reveal 
that the disk evolution is more complicated than the age 
profiles indicate. In all cases, there are few truly young stars 
at large radii; beyond 2 kpc, many fewer stars form after 5 
Gyr. Thus, in contrast to the classical "inside-out" formation 
expected for massive disks, the characteristic radius for star 
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Table 1. Simulation data 





baryon 


V200 










^hreak 




(Mo) 


fraction 


(km s^-'-) 


(km s~^) 


(Mq) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


109 


10% 


15 


7.4 


2.34 X 10"^ 


0.37 


0.30 


0.71 


0.80 


2.5 X 10^ 


10% 


20 


11.4 


2.15 X 10^ 


0.67 


0.53 


0.90 


1.24 


5 X 10^ 


10% 


25 


15.1 


7.86 X 10^ 


0.97 


0.63 


0.54 


1.33 


8.6 X 10^ 


10% 


30 


20.1 


2.20 X 10** 


0.63 


0.84 


1.95 


1.24 


1.4 X 10^0 


0.5% 


35 


15.6 


5.23 X 10*^ 


0.65 


0.47 


2.09 


1.06 


1.4 X lO^o 


1% 


35 


16.8 


1.72 X 10^ 


0.84 


0.61 


2.37 


1.15 


1.4 X IQio 


10% 


35 


29.9 


3.86 X 10** 


0.35 


0.91 


3.45 


1.24 



formation shrinks with time for these smaller systems. Our 
high mass model galaxies do show an early period of "inside- 
out" growth during which the small age gradient seen in 
Figure |4] is imprinted at large radii. However, after 5 Gyr, 
they too experience the same contraction in the star forming 
region of the disk that is seen in lower mass galaxies. 

The contracting star formation results from the interac- 
tion between gas pressure and the star formation density cri- 
terion. The gas disks of these small systems are partially sup- 
ported by gas pre ssure in addition to classical angular mo- 
mentum support llRhee et al.l [20041 : [Valcnzuela et al. 2003; 
iKaufmann et al.ll2007l ')~ Therefore, as gas is consumed in the 
center due to initial high star formation rates, the gas disk 
contracts radially until pressure support is re-established. 
The resulting decline in central star formation rate may also 
reduce the turbulent velocity of the gas, further reducing the 
pressure support. This process leads to radial shrinking of 
the gas disk such that high gas densities are maintained in 
the center, at the expense of the outer regions. The resulting 
evolution leads to star formation that is systematically less 
radially extended as time goes on. 

We can verify that the age gradient observed in our 
simulations is established at formation rather than through 
stellar migration. Figures [7] and [8] show the radius where 
stars formed as a function of time. The overall shape of 
Figures [7] and |8] closely resembles Figures [5] and |6l both 
sets of figures show shrinking of the radial envelope within 
which stars form. Formation location is therefore a primary 
factor in determining the final age gradient, rather than the 
degree of stellar migration. One aspect of Figures [5] and [6] 
that differs greatly from Figures [7] and [5] is that many stars 
finish the simulation beyond where they form. This outward 
movement is investigated in more detail in the next section. 

Each of the models shows two distinct spatial regions. 
At small radii, the star formation is abundant and nearly 
continuous. At large radii, there is less star formation and 
what there is is episodic. During these episodes of elevated 
star formation, gas temporarily reaches star forming densi- 
ties in the outer region of the disk, typically in spiral pat- 
terns. Once the stars have formed, the supernova feedback 
disperses the dense gas and shuts down star formation. 

In models with lower stellar mass, the episodic bursts 
are scattered at regular time intervals spaced 300 Myr 
apart. This timesc ale is similar t o the star formation in- 
terval reported in IStinson et al.l (|2007l ). and is related to 
the free fall time at the center of halos. Higher mass galax- 
ies also show similar episodic star formation, but only out- 
side the inner stable star forming disk. At large radii, the 
episodic star formation has a longer timescale of ~ 1 Gyr 



due to the lower characteristic densities associated with 
longer free fall times. This episodic star formation can be 
explained with a delay different ial equation as shown by 
lOuillen fc Bland-HawthornI ([2008). 

For a more direct comparison with observations. Fig- 
ure]^ plots model CMDs showing the expected stellar popu- 
lations in several radi al bins. To cr eate the CMDs, we used 
StarFISH version 1.1 (IHarris fc Zaritskv 200l]) to populate 
a set of isochrones l|Girardi et al.l |2002| ) based on a user- 
supplied star-formation history (SFH) set by the mass, age, 
metallicity (-2.3 < [Fe/H] <-l after 13.5 Gyr), and position 
of the simulated star particles. To mimic observing condi- 
tions, we adopte d artificial star tests fr om the ACS Nearby 
Galaxy Survey l|Dalcanton et al] |2008| ) assuming that the 
simulated galaxy is located at a distance of 500 kpc. 

The CMDs of the simulations show a prominent young 
main sequence in the central 2 kpc that disappears at large 
radii. The absence of the young main sequence outside 2 
kpc agrees with the age gradient shown in Figure [J] Figure 
[9] also shows that any main sequence stars in the halo are 
more than 1 magnitude fainter than the red clump and hor- 
izontal branch, such that only the deepest photometry in 
the nearest galaxies would show evidence of this intermedi- 
ate age population. The age gradient apparent in Figure [9] 
is consistent with the deep CMDs presented in Figure 2 of 
Hidalgo et al. (2003b) suggesting that isolated models can 
produce realistic stellar structure. 



4.3 Stellar Migration 

While the location of star formation largely determines the 
age gradient and stellar structure of the galaxies, there are 
some features that are due to the motion of stars after they 
form. Figures [To] and [TT] highlight these differences by show- 
ing how far the stars have migrated radially. 

The mean migration (dashed line) in the constant 
/6=0.1 models shown in Figure [TO] shows little deviation 
from zero (solid line), although many individual particles 
move a large distance. For the lowest mass halos, most of 
this movement is in the outward direction and is concurrent 
with star formation episodes. 

To explore this behavior. Figure [12] shows the stellar 
apocenters as a function of their velocity at the time of for- 
mation. Compared to the vertical line at zero velocity, stars 
to the right formed from gas that was moving outwards, 
while stars to the left formed as the gas was collapsing. 

There is a noticeable trend for stars with larger apocen- 
ters to form with larger outward velocities, particularly for 
the lowest mass 15 km s"'^ halo. Physically, these halo stars 
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Figure 5. The final distribution of where stars are versus their formation time for the four models with /t,=10% and varying masses 
represented by their V20Q velocities. The plot replaces clusters of dots with contours in particularly dense regions. The contours are are 
based on a density map with cells that are 0.3 kpc by 0.1 Gyr. The contours are scaled between 10 levels from a minimum of 20 stars 
per cell {black) to the maximum of stars in a given cell {white inside black). Younger ages are to the right. 
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Figure 6. The same as Figure [5] except for the models with 11200 = 35 km s ^ and /{, = 0.005, 0.01, and 0.1. 
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Figure 7. The distribution of where stars are formed versus when they formed in the four models with /(,=10% and varying masses 
represented by their D200 velocities. The contours are are based on a density map with grids that are 0.3 kpc by 0.1 Gyr. The contours 
are scaled between 10 linearly spaced levels from a minimum of 20 stars per cell to the maximum of stars in a given cell. 
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Figure 8. The same as Figure [7] except for the models with ?J200 = 35 km s ^ and = 0.005, 0.01, and 0.1. 



J 



1 tl kpc 
I I I I I I I ' I I M I I I ' I II I I I ' 




Feedback and the Formation of Dwarf Halos 9 

1 < T < 2 kpc T > 2 kpc 



1 " 1 1 1 1 1 1 1 ' 1 1 1 


1 1 1 ' |i 1 1 1 I ' 1 1 1 


■ - 


J 




S 






■■ 














-JJ] -03 OuCi DJ JjQ ]J 2.D 



a" 



1 1 1 1 1 1 1 1 



1 1 II 1 1 1 ' 




-IJ] -Cl3 OuCi DJ JjCi ]J 2.D 
Vi.-[.„ 



Figure 9. Synthetic color magnitude diagrams of tlie 15 km s~^, /i,=0.1 simulated dwarf galaxy, shown in three radial bins. Magnitudes 
are in the HST/ACS filter system. The contour colors, from light to dark, represent point densities of 30, 50, and 70 points per 0.1 X 0.1 
mag bin. 



formed from outward flowing gas that was accelerated by 
supernovae blastwaves and shocks against infalling gas. The 
shock produces high densities and triggers star formation. 
These stars are thus launched on radial orbits that create 
an extended stellar halo and a positively skewed velocity dis- 
tribution. Conversely, few stars form during the infall phase 
during which shocks are weaker or absent. No stars are ob- 
served to be created with sufficient energy to become un- 
bound from the dwarf galaxy. Once the dwarf interacts with 
other satellites, some of these stars may become unbound. 

Figure [T^] suggests that shocks and instabilities in the 
SN-driven wind may be necessary for the formation of the 
most extended halo stars seen in dwarfs. To verify that the 
large radius, high initial velocity stars are indeed due to 
feedback, we have rerun the 15 km s~^ simulation with feed- 
back turned off. In this simulation, no extended halo forms. 
Instead, the distribution of apocenters is centered at zero ra- 
dial velocity, consisting of disk stars formed from gas moving 
in stable circular orbits. 

The signature of halo formation due to supernova feed- 
back is much reduced in higher mass galaxies, which form 
long-lived stable stellar disks. In these cases, disk stars form 
on circular orbits and appear near the zero initial radial ve- 
locity line in Figure [T^ T hese stars may mig rate outwards 
through disk instabilities l|Roskar et aLllioOa ). but are un- 
likely to swell into a 3-dimensional halo. 

Figure [T5] also suggests that in the lowest mass galaxies, 
one expects a kinematic signature where the outermost halo 
stars are all preferentially on radial orbits. However, these 
stars can be strongly influenced by tidal effects, possibly 
making this signature difficult to detect observationally. 



4-3.1 Stellar Disk Meandering 

Unlike in the /i,=10% models, the low ft, models show an 
offset between the mean initial and final radii as shown in 



Figure [TT] This offset develops because while the halo poten- 
tial well is deep enough that gas cools onto a disk, the gas 
is pressure supported and not dense enough to form stars 
outside the central region. Thus, in these low models, the 
gas disk is significantly more massive than the stellar disk 
and dominates the disk dynamics, as shown in Figure 1131 
The stellar disk then meanders in response to the dominant 
influence of the gas. As the stellar disk meanders. Figure [TT] 
shows that it scatters stars into the halo. Stars that formed 
at the beginning of the simulation are even further from the 
center than those that formed at the end of the simulation. 
The meandering disk therefore provides another mechanism 
for reinforcing the age gradient imprinted when the stars 
formed. A similar mechanism is apparent in the cosmol ogi- 
cal dwarf galaxy simulation of iMashchenko et al.l (|2008l ) . 

Figure [T?] shows the effect the meandering stellar disk 
has on the age proflle. It compares the final age profiles from 
Figure [4] with age profiles as a function of formation radius. 
Figure [2] thus shows the age gradient that would have been 
observed in the absence of stellar migration compared to 
the final age profile. For the /t = 10% simulations, the stars 
do not meander and their age profile is determined almost 
entirely by where they form. The 35 km s~^ models with 
/(,=0.5% and fb = l% show declining mean stellar age at large 
radii at formation (dashed line), but constantly rising mean 
age in the present day age profile (solid line). The sloshing of 
the disk is therefore capable of changing the sign of the age 
gradient at large radii, particularly when the baryon fraction 
is low. 



4-3.2 Low Temperature Cooling 

Because the disk sloshing is due primarily to high gas mass 
fractions, one wonders if such gas rich disks would be present 
if we had not neglected low temperature cooling (below 10* 
K) . Our neglect of molecular and metal line cooling leads to 
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Figure 10. The distribution of liow much stars moved versus when they formed in the four models with /i,=10% and varying masses 
represented by their D200 velocities. The contours are are based on a density map with grids that are 0.3 kpc by 0.1 Gyr. The contours 
are scaled between 10 linearly spaced levels from a minimum of 20 stars per cell to the maximum of stars in a given cell. The dashed 
line represents a lack of movement from the stars formation location. The dotted line represents the mean movement of particles at the 
various formation times. 
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Figure 11. The same as Figure [TOl except for the models with V200 = 35 km s ^ and varying /(,. 
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Figure 12. The formation radial velocity of stars compared with their apocenter (maximum radial extent) for three representative 
model galaxies. Contours are linearly scaled to the density of star particles in the 1 km s~^x 0.5 kpc size bins. A dashed line is drawn 
at zero velocity to demarcate outward going forming stars from inward moving stars. The mean formation velocity for each apocenter is 
plotted as the gray line. The left panel shows the 15 km s~^ model with only 10^ gas particles. The center panel shows what happens 
when the same 10^ particle 35 km s~^ model is run without supernova feedback. The right panel shows the plot for the 35 km s~^, 
/j,=0.5% model run with 10^ particles. We note that the 15 km s~^ model run with feedback shows the largest fraction of stars formed 
with positive outward velocities (to the right of the dotted line). The 35 km s~^ model shows slightly more stars forming with outward 
velocities than inward. The supernova feedback is a constant energy and it is able to move stars the farthest out of the shallowest 
potential wells. 
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Figure 13. The star to gas mass ratio in the stellar disk. The 
"stellar disk" is defined using the star with the largest R and z 
from the center of mass of the star particles. This definition of the 
disk is conservative as the gas disk stretches beyond the stellar 
disk in our simulations. Even limiting the definition of the disk as 
we do, the larger gas mass in the /i,=0.5% and /i,=l% models is 
apparent. This larger mass explains why the stellar disk meanders 
through the gas disk. 
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Figure 14. The mean stellar age versus formation radius. This 
shows the age profile if the stars did not move at all after they 
formed. The mean is taken of stars grouped together into 40 radial 
bins with equal numbers of particles. 



an artificial temperature floor that supports a large l o w den - 
sity gaseous disk, as explained in iKaufmann et al] (|2007h . 
However, low mass galaxies contain few metals, so low tem- 
perature cooling should be suppressed making our neglect 
of these processes less important than for massive galax- 
ies. Low metallicity also suppresses molecular cooling, for 
two reasons. First, the formation of molecules is inefficient 
when dust grains are absent, greatly reducing the mass of 
molecular gas when the metallicity is low. Second, molecular 



cooling is dominated by emission lines from metal-bearing 
molecules, leading to more than a factor of ten reduction in 
the coohng rates when metals are absent. Moreover, there is 
clear evidence that even with atomic cooling alone, our sim- 
ulations contain gas that is cold and dense enough to form 
stars. Observationally, extended HI disks are a normal fea- 
ture of dwarf galaxies (jSwaters et al.ll2002l ). Combined, this 
suggests that our neglect of low temperature cooling may 
not be a significant limitation in these simulations. 
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4.4 Broken Density Profiles 

The density profile of the high resolution 15 km s^^ model in 
Figure |3]shows a surprising amount of structure for a galaxy 
that formed in isolation without any hierarchical buildup of 
material. We can relate the structure to the different regions 
of star formation. The innermost volume (r < 1 kpc) cor- 
responds to the radius inside which stars continually form. 
The next region (1 < r < 2 kpc) corresponds to the volume 
inside which heavy star formation occurred only in the past. 
Beyond 2 kpc, there is a modest increase in stellar surface 
density slope. This corresponds to the region in which no 
star formation happened in situ, but into which stars were 
ejected. Using these different regimes of star formation, it is 
possible to see how a dwarf galaxy that forms in isolation 
can create a structured stellar distribution. 

4.5 Halo Shapes 

The shapes of the stellar halos that formed in our simu- 
lations were influenced by how they were produced. A halo 
dominated by stars that formed in supernova shocks is likely 
to be more spherical, whereas one that formed from disk 
shrinking or sloshing is likely to be more flattened like a 
disk. 

Figure [12] shows the axial ratio of the outermost 13% 
of stars as a function of stellar mass (left panel) and of halo 
mass (right panel). There is a general trend for systems with 
low stellar masses and/or halo masses to have more spher- 
ical distributions. The trend with halo mass is due to the 
increasing importance of supernova shock-driven halos in 
galaxies with shallower potential wells. In the more massive 
halos, the formation of the old stellar halo is dominated by 
processes that operate in disks, leading to flatter axial ratios. 

At fixed halo mass, systems with larger baryon fractions 
have fiatter stellar halos because their disks are more ex- 
tended. We note, however, that even the most massive high 
baryon fraction galaxies with the flattest halo axial ratios 
are still rather round in their outskirts. At the mass scales 
explored in these simulations, disks exert significant pres- 
sure support, which leads to thicker star fo rming disks than 
those typical in higher mass galaxies (e.g. iKaufmann et al] 
|2003)). 



5 CONCLUSIONS 

Recent, deep observations of regions surrounding Local 
Group Dwarfs have shown structure in their stellar popu- 
lations in the form of old stellar populations at large radii 
and broken exponential surface density profiles. In larger 
galaxies, these extended stellar halos are thought to be the 
result of hierarchical merging. In our simulations, we find 
that even isolated galaxies are able to form extended halos, 
explaining why so many dwarfs have halos, even when the 
expected masses of the merging sub-halos are thought to be 
too low to form stars. 

In our simulations, three mechanisms contributed to an 
age gradient in the outer galaxy. First, the envelope inside 
which star formation occurred contracted as the supply of 
gas diminished and pressure support was reduced. Second, 
gas-driven disk instabilities caused bulk motion of the stellar 
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Figure 15. Halo shape as a function of stellar mass (left panel) 
and halo mass (right panel) after 13.5 Gyr of evolution of the 
models. The halo shape c/a values come fro m the outermost 
13% of th e stars . The v alue of c/a is defined in iKat j lll99ll ) and 
[Pcbattist a et al as the ratio of moment of inertia between 

the most minor axis and the major axis. The size of the circles 
corresponds to the baryon fraction (/;,) of the model. The stars 
on the outskirts of the smallest galaxies have the least flat shape. 



disk over time so that stars that formed the longest ago 
are the farthest from the current center. Third, some stars 
were ejected far out into the halo after forming in supernova 
driven shocks. The ejection was somewhat more effective at 
early times because of the active early star formation. 

The importance of the three different mechanisms var- 
ied with mass and baryon fraction of the model. For exam- 
ple, the contracting star formation envelope occurred in all 
of our models, but began earliest in the lowest mass models. 
The higher mass models were able to initially grow disks 
from the "inside-out" using their large gas reservoir. How- 
ever, dwindling gas supplies halted the steady disk growth 
and star formation contracted back towards the galaxy cen- 
ters, in a similar manner to the lower mass models. 

Bulk motions of the stellar disk were the most notable 
in the models with low baryon fraction, which replicate the 
effect of reionization by removing baryons from low mass 
halos. Lowering the baryon fraction leads to low density 
gas disks that can only form stars at their center. The low 
baryon fraction disks therefore remain gas rich and are more 
prone to instabilities that lead to sloshing of the embedded 
stellar disk. 

The ejection of stars formed in supernova-driven shocks 
is most prominent in the lowest mass models. At these mass 
scales, the depth of the potential well is shallow enough 
that supernovae can drive significant outflows, making stel- 
lar ejection more likely. 

Each of these three scenarios played a role in producing 
a stellar distribution with a remarkable amount of struc- 
ture. Although we have not ruled out the possibility that 
the structure in dwarf stellar halos is the result of hierarchi- 
cal merging, tidal interactions, or cold accretion along fila- 
ments, it is possible to produce extended stellar structures 
without merging. 
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